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Dynamics and (de) localization in a one-dimensional tight-binding chain 
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A simple tight-binding model is used to illustrate how the time dependence of a state vector 
can be obtained from all the eigenvalues and eigenvectors of the Hamiltonian. The behavior of the 
eigenvalues and eigenvectors is studied for various parameters and allows us to study scattering-like 
events, impurity states, and localization in disordered systems. 



I. INTRODUCTION 

For most physics students the Schrodinger representa- 
tion of quantum mechanics is more appealing than the 
matrix or Heisenberg representation. The reason is due 
in part to the fact that the Schrodinger representation is 
more suited for visualization.^ Also most students are 
more accustomed to thinking in terms of functions rather 
than in terms of abstract eigenvectors. 

In this article I discuss a simple but versatile tight- 
binding Hamiltonian whose eigenvalues, eigenvectors and 
dynamics can be obtained and easily visualized in the 
matrix representation. The solution of the physical sys- 
tem represented by the Hamiltonian is obtained numer- 
ically by using a short program that is given in Ap- 
pendix A. The code allows for visualization of the spec- 
trum of eigenvalues and eigenvectors and also the dynam- 
ics generated by the Hamiltonian. 

The model Hamiltonian of a one-dimensional chain 
with nearest neighbor couplings is discussed in Sec. ^ 
and a complete set of states describing the single parti- 
cle dynamics is introduced. Section HTT1 briefly describes 
how the Hamiltonian generates the time dependence of 
an initial state vector and demonstrates that the eigen- 
values and eigenvectors define the dynamical behavior 
of the system. In Sec. IIVI I discuss several applications 
that are of pedagogical interest. In particular, the dy- 
namics of a state initially localized on a particular site is 
studied. The dynamics is correlated with the spectrum of 
the Hamiltonian and properties of its eigenvectors (local- 
ized versus delocalized). Several aspects of the problem 
that are of research interest are discussed, including the 
problems of localization and conductance in chain-like 
molecules, such as DNA. 



II. DEFINITION OF THE PROBLEM 

Consider the Hamiltonian: 

H = Eic\ci + ^ T i \4+i c i + c l c i+i] > (!) 

i i 

where cj and Ci are the creation and destruction opera- 
tors of a particle on site i, respectively. The Hamiltonian 
represents a chain of sites denoted by indices i; a particle 



can hop from one site to another due to the nonvanish- 
ing values of Ti which are often called hopping matrix 
elements. 

We restrict our attention to a single particle (or, more 
generally, an excitation) propagating through the chain. 
The dynamics can be described in a position-occupation 
basis, that a basis of states denoted by {\i),i — 1, . . . , N} 
such that 



\i) = 00... !,...!! . 



(2) 



where i denotes a particular site occupied by an excita- 
tion and N denotes the total number of sites. The action 
of c-operators in this basis is simple: 



4|0) = \i) 

= |0) 

Ci\j) = 0, i j^j, 



(3a) 
(3b) 
(3c) 



where we consider only singly occupied states; the vac- 
uum is denoted by |0). In the basis that is restricted 
to singly occupied states, we can equivalently (and with- 
out reference to creation/destruction operators) repre- 
sent the Hamiltonian as 

h = J2 m (<i + E T « [n + !) w + n> + !|] ■ w 



The position-occupation basis does not diagonalize the 
Hamiltonian in Eq. JIJ, except in the trivial case of Tj = 
for all i. However, it is the basis that is simplest concep- 
tually and most easy to visualize. 

Equation is a simplified tight-binding Hamiltonian 
and is discussed in many textbooks on condensed matter 
physics (see for example, Rcf. 3). Its matrix representa- 
tion is easy to construct. The simplest and numerically 
most feasible way is to consider a tridiagonal matrix in 
the position-occupation basis (Eqs. an d ©) as 
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The Hamiltonian is represented in a basis of singly occu- 
pied states; that is, we consider only this subset of states 
of Fock space (which includes multiply occupied states; 
for example, we might want to consider the dynamics 
of two excitations). Note that the Hamiltonian defined 
in Eq. Q cannot induce transitions between the Fock 
subspaces corresponding to a different total number of 
excitations. 

Periodic boundary conditions are not imposed in the 
matrix representation in Eq. JSJ. Periodic boundary con- 
ditions would require nonvanishing (1, AT) and (N, 1) el- 
ements, that is, the upper right and the lower left corner 
of the matrix. The matrix in Eq. JSJ) can be easily set 
up and diagonalized numerically, ft is enough to specify 
only two arrays of real numbers, one of length N which 
contains the diagonal values of the matrix, and the other 
of length N — 1 which contains the elements of the Hamil- 
tonian matrix along its first subdiagonal. 

III. TIME DEPENDENCE OF STATE VECTORS 

Let us assume that all the eigenvalues tu and eigen- 
vectors of the Hamiltonian H are known (if |£fc) = 
We further assume that the system is at time 
t = 0, in some known or prepared state \^(t = 0)). The 
state \fy(t = 0)) can be projected onto the basis of eigen- 
vectors of the full Hamiltonian: 

= 0)> = 5>*l&>, (6) 

where the projection coefficients a>k are given by 

a fc = <&|*(i = (>)>, (7) 
because the |£/c) are assumed to be orthonormal, that is, 

(£m|£n) = 5m,n- (8) 

Let us denote the time evolution operator by U (t) ; that 
is, U(t) acts on an arbitrary state \^(t — 0)) and evolves 
it to the state | 

um(t=o)) = m)). (9) 

The action of the time evolution operator on the eigen- 
vectors of the problem is trivial^: 

U(t)\{ k ) = exp(-iekt/fiMk), (10) 
which implies that 

l*(*)) = 53 afeex PHefe*A)l&) (H) 

k 

= X)«fcl*(< = 0)>exp(-te fc t/R)|a). (12) 

k 

We assume that the wave function (or state vector) is 
initially chosen to be localized on a particular site I of the 



chain, that is, \ty(t = 0)) = |Z). We ask about the prob- 
ability that after some time t the excitation is on some 
other site to. To obtain this information, the state vec- 
tor ^(t)) must be projected onto the position-occupation 
basis, that is, we should calculate N projections 6 m ,z(i) 

b m ,i{t) = (m|*(t)) = $3<&|i) H&>e*p(-*e*W 

k 

. (13) 

The quantity |& m ,z(i)| 2 is the probability that at time t, 
the excitation initially localized on site I is found on site 
to. A pedagogical account of the definition of probability 
current in tight-binding problems can be found in Ref. la. 



IV. APPLICATIONS OF THE MODEL 

As mentioned, the matrix in Eq. JSJ can be diagonal- 
ized numerically. One possible way of doing so is de- 
scribed in Appendix A, which lists a sample program to 
set up the Hamiltonian matrix and diagonalize it. The 
result of this numerical procedure is an array of N Hamil- 
tonian eigenvalues e k and a, N x N matrix (or 2D array) 
of eigenvectors which provides a complete solution 
of the problem, including its time dependence. In the 
following I present several applications of the code. 



A. Regular chain, propagation of the initially 
localized state 

I first consider a completely regular chain, that is, a 
chain in which Ei = E and Tj = T for all i. The char- 
acteristic energies (E) and the hopping matrix elements 
(T) are set to E = e and T = — O.le, where e denotes 
the energy scale. Note that the chosen energy scale also 
fixes the characteristic time scale, which is given by %/e 
(see Eq. 11131 ). The results of this calculation are given 
in Figs, ^and |21 It is clearly observed how the initially 
localized state delocalizes over many sites and eventu- 
ally hits the right end of the chain, reflecting from it. 
The horizontal axis in the plots is the site index which is 
in principle unrelated to any characteristic length. The 
spatial dimensions are hidden in the hopping (or over- 
lap) matrix elements Ti; the more separated the atomic 
orbitals Ei are, the smaller their overlap and the corre- 
sponding hopping matrix element^ 

Figure [21 displays the eigenvalues of the Hamiltonian 
and three characteristic eigenvectors. Note the fast oscil- 
latory behavior of the eigenvectors for high-energy states 
(insets in (c) and (d) panels of Fig. 0). The eigenvalues of 
the Hamiltonian are indistinguishable from the analytical 
solution, which is given by& 

e k =E + 2Tcos[kir/(N + l)}, (14) 

where k is the eigenvalue index (k = 1, . . . , N), that is, 
a cosine band of states of width W w 4T = 0.4e. The 
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FIG. 1: Time dependence of the probability distribution 
|&m,i| 2 for a completely regular chain. The excitation was 
initially localized on the first site, that is, \9(t = 0)} = |f). 
The time is measured in units of Ti/e and the number of sites 
N = 500. 

exact solution for the periodic tight-binding chain ism 

e(k) = E + 2Tcos(2kn/N), (k=l,...,N) (15) 

Note the extra factor of 2 in the argument of the cosine 
compared to Eq. (|14|l and that the eigenvalue e& has been 
rewritten as e(fc), so that it appears as a function of the 
eigenvalue index. In the periodic case, it makes sense to 
characterize the eigenvalues by the wavevector, that is, k 
becomes more than an eigenvalue index and has a direct 
interpretation in terms of the characteristic wavelength 
for each eigenstate. This issue is discussed in more detail 
in Appendix IbI 

Note that Eqs. I|14|l and (|15|l are not very different 
for large N. In the limit of infinite N, the band width, 
defined as the difference between the largest and smallest 
eigenvalues, is the same in both cases, as well as the 
density of states defined as 

p(e)=J2S(e-e k ). (16) 
ft 

There is one important difference, that is, the double 
degeneracy of states given by Eq. ifTBTl for states with 
k% = I and &2 = N — I, I — 1, . . . , N — 1, which is not the 
case in Eq. (|THl . 



FIG. 2: Eigenvalues and three characteristic eigenvectors of 
a regular chain, (a) Hamiltonian eigenvalues, (b) eigenvector 
of the lowest energy eigenstate (e^ = 0.8e), (c) eigenvector 
of eigenstate with energy = 0.99937e (250th eigenstate), 
and (d) eigenvector of eigenstate with energy = 1.2e (high- 
est energy eigenstate). The insets in the (c) and (d) show 
the eigenvector magnitude between the 100th and 120th sites 
(these are blowups of the eigenvector magnitude in the regions 
bounded by the two thin rectangles in the main plots). 



B. One defect link in a chain, simulation of 
scattering 

In this subsection, a special link is introduced between 
the 250th and 251st sites in a chain with total of N = 500 
sites so that T250 = — 0.2e, and all other links are the 
same as before, Ti = — O.le for all i, i ^ 250 (the orbital 
energies are Ei = e for all i, the same as in previous 
subsection). This modification of the hopping matrix 
will allow us to study the effects of the impurity link on 
the eigenvalue spectrum and propagation of the initially 
localized state. 

The evolution of a state vector initially localized on the 
first chain site is displayed in Fig. [3J Note how part of 
the probability density is reflected from site 250 and 251, 
while the other part continues its propagation toward the 
end of the chain. 

The eigenvalue spectrum is shown in Fig. 0] Note the 
appearance of two states that detach from the band. The 
eigenvectors of two special states are also displayed in 
Fig. The two states separated from the band corre- 
spond to the excitations that are localized on the special 
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FIG. 3: Time dependence of the probability distribution for 
a regular chain with a special link inserted between the 250th 
and the 251st site of the chain (denoted by dashed lines in 
the plots). The time is measured in %/e. 



sites 250 and 251, that is, these states are related to ex- 
citations of impurity sites. 



C. Anderson's diagonal disorder Hamiltonian 

We now study a case in which randomness is intro- 
duced in the Hamiltonian matrix. The orbital energies 
Ei are given random values in a band of width W. Be- 
cause these numbers are along the diagonal of the Hamil- 
tonian, the model is said to have diagonal disorder. The 
sub-diagonal matrix elements are the same as before 
(Tj = T = — O.le). There are a number of interesting 
issues related to this model, one of which is called An- 
derson localization. The suitably modified code in Ap- 
pendix A can be used to study the Anderson's diagonal 
disorder Hamiltonian. The number of sites is increased 
to N = 900 and the orbital energies along the diagonal 
are given uniform random values in the interval [e, 2e], 
that is, W = e. 

One of the features of this problem is that all the eigen- 
states of the Hamiltonian are localizedAS The eigenvalue 
spectrum and three characteristic eigenvectors are shown 
in Fig. [5] The fact that all the eigenvectors are localized 
has a profound influence on the propagation of an ini- 
tially localized state. Figure displays an evolution of a 



FIG. 4: Eigenvalues and the two eigenvectors corresponding 
to impurity excitation states, (a) Hamiltonian eigenvalues, 
(b) eigenvector of the lowest energy eigenstate (e^ = 0.75e), 
and (c) eigenvector of highest energy (e& = 1.25e). The two 
sites that are connected with an impurity link are denoted by 
dashed lines in (b) and (c). 



state vector initially localized on the first chain site. Oc- 
cupation probabilities for the first 30 sites of the chain 
arc presented; the site occupation probabilities are negli- 
gible beyond site 10. Note that the times shown are very 
long (t = 10 w h/e). We conclude that the propagation of 
the excitation through the chain is not only slow, but is 
effectively blocked - the excitation remains localized in 
the vicinity of the site at which it was initially created£& 
The blockage of the excitation propagation is related to 
the fact that the projection coefficients (afc in Eq. ©) of 
the initially localized state have a significant magnitude 
only for several eigenstates whose localization on the first 
site is nonvanishing. One of those eigenvectors (the one 
whose maximum magnitude is on the first site) is very 
similar to the initial state vector, and its projection coef- 
ficient is the largest and close to 1. Thus, the initial state 
is almost an eigenstate and its evolution is thus slow. Be- 
cause the projection coefficients at on eigenvectors that 
are localized on sites that are very distant from the first 
chain site are close to zero, the propagation through the 
chain is essentially blocked. The total number of the sites 
that become occupied during the evolution (about 5 to 10 
as is seen in Fig. |SJ) is related to the typical localization 
width of the eigenvectors. 
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FIG. 5: Eigenenergies and three characteristic eigenvectors of 
the Anderson's diagonal disorder Hamiltonian. 



D. Recent research on chains in the tight-binding 
scheme 



FIG. 6: Time dependence of the probability distribution for 
a chain with random orbital energies chosen from the interval 
[e, 2e]. The time is measured in h/e. Note that only the first 
30 sites of the chain with 900 sites are represented, because the 
site occupation probabilities are negligible beyond the 10th 
site. 



It is interesting to see how the possible correlations 
in the diagonal disorder distribution influence the na- 
ture of the eigenvectors and th e sp ectrum of eigenvalues 
(see for example, Refs. and llOf) . These effects could 
be studied with suitable changes of the code. The lo- 
calization of electronic states due to disorder may drive 
a metal-insulator transition (or Anderson transition) in 
the system, and one-dimensional models similar to the 
one discussed in this article are useful in the study of 
random binary alloys^i In this case there are only two 
characteristic orbital energies and their appearance in the 
chain is random (for example, the appearance of the A 
orbital occurs with probability p, while the appearance 
of the B orbital occurs with probability 1 — p). For a 
random binary alloy, the system again exhibit the local- 
ization of eigenvectors, but correlations in the disorder 
may introduce resonant states for which there is perfect 
electron transmission through the system. 

The electrical conduction of biological polymers, DNA 
in particular, has also attracted much attention>Aiii2 An 
electronic coupling induced through the overlap of tt z or- 
bitals perpendicular to the planes of the stacked base 
pairs in double-stranded DNA can be simulated using the 
simplified one-dimensional model in Eq. Q , although the 
realistic situation is much more complicated due to the 



influence of vibrations on the distances between the base 
pairs and the importance of the electronic structure of 
the DNA backbonei 12 i 13 i 14 i 15 Models similar to the one 
studied in this article are used to study the conduction 
properties of DNA moleculesiii 

The model and the program in Appendix ^ can be 
easily modified to study these problems. Correlations 
in the diagonal disorder can be studied, as well as the 
introduction of disorder along the sub-diagonals (hopping 
matrix elements). Localization effects can be studied in 
cases when the width of the diagonal energy band W 
is much smaller than the value of the hopping matrix 
elements T. The opposite limit was studied in Sec. lIVCl 



APPENDIX A: A PROGRAM FOR THE STUDY 
OF THE PROBLEM 

The following code describes the numerical solution to 
the problem. Due to the brevity of the code, it is listed 
here along with comments. The source code can be com- 
piled as g77 -o chainl chainl.for -llapack or f77 
-o chainl chainl.for -llapack 

depending on the Fortran compiler installed. Note that 
the code should be linked to the LAPACK library of 
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routines^*. 

1 program tbald 

2 double precision E(5000) , T(5000) , XI(5000,5000) , 
& flag, work(9998), thresh, time, ReB, ImB 

3 open(l, file=' eigenvalues. dat' , st atus=' unknown ' ) 

4 open(2, file=' eigenvectors .dat ' , st atus =) unknown ' ) 

5 N = 500 

* Number of chain sites 

6 thresh = 0.0D0 

* Threshold for random function 

7 do i=l, N 

8 E(i) = 1.0 

* E - diagonal of hamiltonian matrix - orbital 

9 flag = randO 

10 T(i) = -0.1 

* T - subdiagonal of hamiltonian matrix - hopping 

11 if (flag .It. thresh) then 

12 E(i) = 1.0 + 1.0*rand() 

13 endif 

14 enddo 

15 call DSTEVC'V ,N,E,T, XI, 5000, work, info) 

16 do i=l, N 

17 write(l,*) i, ' E(i) 

* On exit from DSTEV, E contains eigenvalues 

18 write(2,*) 'Eigenvector: ', i 

19 do j=l, N 

20 write (2,*) XI(j,i) 

* On exit from DSTEV, XI contains eigenvectors 

21 enddo 

22 enddo 

23 close(l) 

24 close(2) 

25 1 = 1 

* Index of initially occupied site 

26 print *,' Input : time t' 

27 read *, time 

28 open(4, f ile=' timedat . dat ' , status= ' unknown ' ) 

29 do j=l, N 

* j counts the chain sites 

30 ReB = 0.0D0 

31 ImB = 0.0D0 

* Real and imaginary parts 

32 do i=l, N 

* i counts the eigenstates 

33 ReB = ReB + XI (1 , i) *dcos (E(i) *time) *XI ( j , i) 

34 ImB = ImB - XI (1 , i) *dsin(E(i) *time) *XI ( j , i) 

35 enddo 

36 write(4,*) j, ' ', ReB**2 + ImB**2 

37 enddo 

38 close (4) 

39 end 



The total number of sites in a chain (iV) is defined 
in line 5. The Hamiltonian matrix is defined between 
lines 7 and 14. The parameter thresh allows for the 
introduction of random orbital matrix elements (lines 11— 
13); for thresh = the orbitals are regular, while for 



thresh > 1 they are totally random within the ranges 
defined by line 12 (between 1 and 2 energy units). 

The Hamiltonian matrix is diagonalized in line 15 us- 
ing the DSTEV LAPACK routine designed for the diag- 
onalization of tridiagonal matrices. 16 After the diagonal- 
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ization, the jth column of the matrix XI (i, j) contains 
the eigenvector corresponding to the jth eigenvalue of 
the Hamiltonian The DSTEV routine sorts the 

eigenvalues in an ascending order. 

The time dependence of the initial state vector, one of 
the vectors from the {\i}} basis in Eq. (J2J, is implemented 
between lines 25 and 37. In particular, the part of the 
program between lines 29 and 37 implements Eq. (|13|l . 
The initial state is specified by the variable 1 which rep- 
resents the index of the occupied chain site at t = 0. 

The output is written in the file timedat.dat which 
can be plotted separately. This output was used to gen- 
erate Figs.nOl and El The CPU time needed to calcu- 
late the probability distribution does not depend on the 
physical time input (variable time) because the program 
does not propagate a solution in the time domain. All 
the calculations needed for the evolution of the initial 
state vector are performed by the diagonalization of the 
Hamiltonian matrix. This output is used to calculate 
the probability distributions for arbitrary times t (see 
Eq. 0). 

To implement the case of regular chain discused in sub- 
section IIV Al all randomness must be eliminated, and 
thus, the thresh variable in line 6 is set to zero. 

To simulate one special link (impurity or defect) in the 
middle of the chain with N = 500 sites (the case studied 
in subsection liVBfl , we introduce the following piece of 
code in between the lines 13 and 14, 

if (i .eq. 250) then 

T(i) = -0.2 
endif 

The rest of the lattice is in a regular state, that is, line 6 
which controls the amount of randomness in the chains 
is still 

thresh = . 0D0 

For the case of totally disordered chain with N = 900 
sites studied in subsection IIV CI lines 5 and 6 of the code 
should be changed to 

5 N = 900 

6 thresh = 1 . 0D0 



APPENDIX B: PERIODIC VERSION OF THE 
REGULAR TIGHT-BINDING CHAIN 

To obtain a chain that is periodic, we have to connect 
the first and the iVth sites by the hopping matrix ele- 
ments. Hence, the Hamiltonian matrix is no longer tridi- 
agonal because it now contains nonvanishing (1,N) and 
(N, 1) elements. However, the Hamiltonian matrix is still 
symmetric. The periodic chain is easily programmed us- 
ing techniques similar to those described in Appendix lAl 
(the code can be obtained from the author) . For the peri- 
odic case the upper (or lower) triangle of the Hamiltonian 
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FIG. 7: Comparison of the eigenvalues and two characteristic 
eigenvectors for a uniform chain with non-periodic (squares) 
and periodic (circles) boundary conditions. The parameters 
are E = e, T = — O.le, and N = 100. (a) Eigenvalues of non- 
periodic and periodic chain, (b) eigenvectors of the lowest 
energy eigenstates, and (c) eigenvectors of the highest energy 
eigenstates. 



matrix and not only its diagonal and main subdiagonal 
has to be stored. 

The results for a regular periodic chain are shown in 
Fig- E|and compared with the case studied in Sec. IIV Al 
The parameters of the calculation are E = e, T = — O.le, 
and ./V = 100. The results are used to illustrate sub- 
tleties discussed at the end of Sec. IIV Al Note that the 
eigenvalue spectra look indistinguishable, but a closer in- 
spection (inset in Fig. Efa)) reveals that the spectrum 
corresponding to the periodic chain is doubly degenerate 
(except for the lowest and highest energy eigenvalue), 
which is not the case for the nonperiodic chain, in agree- 
ment with the discussion in Sec. IIV Al The eigenvectors 
corresponding to the lowest and highest energy eigenstate 
are also different (which is true for all eigenvectors) al- 
though the energies of these states are similar in the two 
cases, see Eqs. I|14|) and (QSJ). For the periodic chain the 
lowest energy eigenvector can be written in the position- 
occupation basis as 



lowest) = -== 1 1,1, 1,1... 1>, (Bl) 

V iv 
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and the highest energy eigenvector in the same basis is 
highest) = -/=!- 1.1. -1,1... 1). (B2) 
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